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MULTITYPE RANDOMIZED REED FROST EPIDEMICS AND 
EPIDEMICS UPON RANDOM GRAPHS 

By Peter Neal 

University of Manchester 

We consider a multitype epidemic model which is a natural exten- 
sion of the randomized Reed-Frost epidemic model. The main result 
is the derivation of an asymptotic Gaussian limit theorem for the 
final size of the epidemic. The method of proof is simpler, and more 
direct, than is used for similar results elsewhere in the epidemics lit- 
erature. In particular, the results are specialized to epidemics upon 
extensions of the Bernoulli random graph. 

1. Introduction. The randomized Reed-Frost epidemic is a very general 
model for homogeneously mixing SIR (susceptible — ► infective — ► removed) 
epidemic models; see [14] and [18]. That is, while infectious, an infective, 
i say, has probability V{ of making an infectious contact with any given 
susceptible member of the population. The random variables {V^} are as- 
sumed to be independent and identically distributed. Furthermore, the only 
transitions in state are: from susceptible to infective, and from infective to 
removed. Since infectious contacts with nonsusceptible individuals have no 
effect, we do not need to make assumptions about the relationship between 
an infective and a nonsusceptible individual. 

It is trivial to show that the generalized stochastic epidemic (see [5]), 
where individuals have independent and identically distributed infectious 
lifetimes, and while infectious make (potential) infectious contacts at the 
points of a homogeneous Poisson point process, is a special case of the ran- 
domized Reed-Frost epidemic. It was shown in [16] that by constructing the 
epidemic and random graph in unison, the generalized stochastic epidemic 
upon a Bernoulli random graph also satisfies the randomized Reed-Frost 
criterion. 
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The asymptotic and exact final size distribution of the randomized Reed- 
Frost epidemic have been derived in [14] and [18], respectively. The aim 
of the current work is to obtain the asymptotic final size distribution of a 
multitype randomized Reed-Frost epidemic. (The exact final size distribution 
is given in [18].) A full description of the multitype randomized Reed-Frost 
epidemic will be given in Section 2. Throughout we shall take the final size of 
the epidemic to denote the total number of initially susceptible individuals 
infected during the course of the epidemic. The analysis of the asymptotic 
final size distribution involves a method of proof which differs markedly from 
that given in [14]. We utilize an embedding argument based upon [19] and 
[21]; see also [6, 8, 9, 10]. In particular, our results generalize those of [6], 
Section 4. Conventionally (see, e.g., [6]), an individual's type is determined 
deterministically. However, we also consider the case where the allocation of 
individuals to type is random. Thus the model considered in [7] is another 
special case. 

The main motivation for the current work was extending the asymptotic 
results of [16] to epidemics upon more general random graphs/population 
networks. Epidemics upon graphs have received considerable attention in 
recent years; see, for example, [3, 4, 16, 17]. Only in [16] is a central limit 
theorem for the final size of the epidemic considered. Thus in Section 7.3, 
we show that epidemics upon various extensions of the Bernoulli random 
graph satisfy the multitype randomized Reed-Frost criterion, and therefore 
we obtain results which are considerably more extensive than those in [16]. 

The paper is structured as follows. In Section 2 the multitype randomized 
Reed-Frost epidemic is defined, and a construction of the epidemic process 
suitable for analysis is given. In Section 3 a weak law of large numbers re- 
sult is obtained for the proportion of the initially susceptible population 
that are ultimately infected by the epidemic. This is followed in Section 4 
by a branching process approximation to the epidemic. This is appropriate 
when the total number of initial infectives is small. The main result is pre- 
sented in Section 5, and is the derivation of an asymptotic Gaussian limit 
for the final size of the epidemic. While the results proved in Sections 3-5 
are "standard" for epidemic processes, the method of proof, especially in 
Section 5, is novel. We present a simpler proof for deriving the asymptotic 
Gaussian limit for the final size distribution. This proof is applicable for a 
wide range of epidemic models. In particular, the approach taken in Sec- 
tion 5 can be applied to the great circle model of [10] to extend the results 
therein. Throughout Sections 3-5, it will be assumed that the allocation of 
individuals to types is deterministic. Therefore in Section 6 we outline the 
(minor) modifications required to consider random allocation of individuals 
to types. Finally in Section 7 the results are specialized to the models of [6] 
and [7], and epidemics upon extensions of the Bernoulli random graph. 
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2. Description and construction of the epidemic process. Consider a 
closed population divided into m types (or groups). Let a= (a±, 0,2, ■ ■ ■ , a m ) 
and N = (N\,N2, • • • , N m ) denote the total number of initial infectives and 
susceptibles, respectively. [We shall take the statement that a group of in- 
dividuals, S say, is of size s = (s\, S2, ■ ■ ■ , s m ), to be interpreted that S con- 
tains Sk individuals of type k (k = 1, 2, . . . , m).] For i = 1,2, . . . , m, let V, = 
(Vi 1, Vi 2, • • • , Vi t7n ) be an m-dimensional random vector taking values on the 
m-dimensional unit cube [0,1], where and 1 denote the m-dimensional 
vectors all of whose components are and 1, respectively. [Throughout this 
paper, all vectors are row vectors, and inequalities between vectors are to 
be interpreted componentwise, with u < v (u < v), if for 1 < i < m, Ui<Vi 
(and at least one of the inequalities is strict). Furthermore, for any vectors 
u, v, w G M m , [u, v] = {w : u < w < v}.] The infectious behaviors of all indi- 
viduals are assumed to be independent. For i = 1,2, ... ,m, assign to each 
infective of type-i an independent infectivity random vector distributed ac- 
cording to Vj. Then for i = 1,2, ...,m, a type-i infective with infectivity 
random vector V (ijl ) = (%i),i, %i),2, ■ ■ , %i),m) has probability of 
making, while infectious, an infectious contact with a given susceptible of 
type k, independent of whether the individual makes any other infectious 
contacts or not. This gives a natural multitype extension of the randomized 
Reed-Frost epidemic described on page 270 of [14]. Clearly, the total number 
of individuals of each type contacted by a given infective can be dependent. 
This is in contrast to [2], where a multitype version of the Martin-L6f [14] 
model is considered. 

For the asymptotic results derived in this paper, it will be convenient to 
consider a sequence of epidemics {E u } as v — > 00. Consider fixed v. We shall 
assume that a u , and Vy (1 < i < m) are dependent upon v. Let N u = 

N{ + N% H h N^, with N v 00 as v 00. For i = 1, 2, . . . , m, let < = 

— ► TTj > 0, and Q = N J — ► Q > as v —> 00. For i = 1, 2, . . . , m, label the 
initial infectives of type i, individuals (i, — (a\ — 1)), (i, — (a" — 2)), . . . , (i,0), 
and label the initial susceptibles of type i, individuals (i, 1), (i, 2), . . . , (£, iVf). 
Assign to the initial infectives of type i, infectivity random vectors V^^, 
2 )' • • • a")' an< ^ ^ or initially susceptible individuals of type i infected 
during the course of the epidemic, assign infectivity random vectors sequen- 
tially from V^ <+1) , V£ j<+2 ), . . . , V( i><+Nr y 

For each v > 1, the epidemic process can be constructed on a generational 
basis by sampling suitably defined stopping times of a simple multivariate 
counting process, X"(-). The embedding of the epidemic processes (and in 
particular, its final size distribution) into the X^(-)-processes is based on 
the approach introduced in [19]; see also [6, 9, 21]. 
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For v > 1 and t > 0, consider the first [tiV'TI] infectives with infectivity 
random vectors, 

( V (l,l)> V (l,2)> • • • ' V (l,[tiiV"7ri])> V (2,l)' V (2,2)' • • • ' V (m,[t m V"7r m ])}> 

where n = diag(7i"i, TT2, ■ ■ ■ ,vr m ). Then for 1 < i < m, and 1 < j < N", let 
Xijft) = 1, if individual is infected when exposed to the first [tiVTI] 
infectives, and x\ = 0, otherwise. Hence for fixed 1 < i < m and t > 0, 
the {Xij(t)} are identically distributed. For i = 1,2, . . . ,m, and t > 0, let 

*f(t) =E^iX&(t), with X"(t) = (Xf(t),^(t),...,X^(t)). Thus X"(t) 
counts the total number of initial susceptibles that would be infected if each 
individual were exposed to the first [tiV'TI] infectives. 

The epidemic process E u can be constructed using the X l/ (-)-processes 
as follows. Let generation denote the initial infectives, and for k > 1, let 
generation k denote those individuals infected by the infectives in genera- 
tion k — 1. For k > 1, let denote the total number of initial suscepti- 
bles infected by the infectives in generations through to k — 1, inclusive. 
Therefore for k > 0, T u k+1 = X"(T£ + CjQ, where T o = °> and for i > °; 
= T^A^n)" 1 . Thus for k > 1, 

T^^x^T^ + Q^n)- 1 . 

Consequently, the above sequence stops at generation where A;* = min{/c : 
T£ +1 = T v k }. Since N v is finite, k* is well defined. Let T£,(= T^(iV<TI)) 
denote the final size of the epidemic. Then T^, = T k « , and satisfies 

(2.i) T^ = x-(T^ + c)(^n)- 1 . 

Finally in Section 6 we shall consider the case where for each v > 1, 
a 1 ' = and = N u tv u are random vectors, with a" and N u constant. 

The above construction of the epidemic process can again be used with the 
assumption that C, u — ► £ and 7r u — 7r as ^ — > oo. 

3. Weak law of large numbers for the final size distribution. In this 
section we shall consider the limiting distribution of as v — > oo. In par- 
ticular, we show in Corollary 3.3 that 

min{|T^|,|T^-r|}-^0 as ^ oo, 

where r is the solution of an appropriate deterministic model. 

We shall begin by deriving r. For 1 < k < m, let EfV^,] = [i v ik , and 
let cov(V^-, V^ k ) = \"j k . We shall require that for 1 <i,k < m, there exists 
< fJ-ik < oo, such that N v ' [i v ik — > /Ujfc and N v \" kk — > as ^ — > oo. 

For f > 1, let 

V" = {V^ !), V|i j2) , . . . , Vj^+jvj,), V( 2jl) , V£j >2 ), • • • , V (m,< n +A^ )} 
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denote the infectivity random vectors of all individuals in the population. 
For v > 1 and 1 < i < m, let r\ (t) = E[x£ 1 (t)] and <(t|V^) = E[x£ 1 (t)|V 1 ']. 
Let Ny = N v iTi (1 < i < m). Then for 1 < i < m, 

rr(t|v")=E[ X r l i(t)|v i ']=i-n n a-^w 

fc=i i=i 

and 

rT(t) = E^[r?(t|V")] = 1 - -rti) ltk ** ] - 

k=l 

For 1 < i < m and t > 0, let rj(t) = 1 — exp(— SfcLi tk^k^ki)- It is trivial to 
show that for 1 <i <m and t > 0, (t) — ► rj(t) as f — > oo. 
For C > 0, let 

T(C)={T>0:r = r(T + C)}, 

where r(t) = (ri(t), T2(t), . . . , r m (t)) (t > 0). Let cr = 1 — r. Then cr = (o"i, 02, 
. . . , a m ) satisfies 

(3.1) cjj = exp^-^7r fc ^ fci (l + Cfc-o-fc)^ (i = 1, 2, . . . , m). 

Equation (3.1) is identical to [20], equation (4) and [6], equation (2.2). There- 
fore denoting the matrix with elements fj,^ by M, it follows from [20], Lemma 
1, that if M T U is irreducible, and £ > 0, then the equation (3.1) for cr, and 
hence, r has a unique solution in [0, 1]. If £ = 0, then r = is a solution 
of (3.1). Let R denote the Perron-Frobenius eigenvalue of Mil. Then for 
£ = 0, there exists r 6 T(0) such that r > 0, if and only if R > 1. We shall 
let r denote the nonzero solution of r = r(r + £) whenever such a solution 
exists (i.e., £ > or R > 1). 

Before analyzing the X^(-)-processes (Lemma 3.2), and hence T^, (Corol- 
lary 3.3), we develop two useful bounds for the covariances of {Xij(')} 
(1 < i < m). The bound (3.3) will not be required until Section 5. 

Lemma 3.1. For l<i<m, 1 < j < I < N? and for all t, u > ; 

m 

(3.2) 0<cov( X r >i (t),xr ) z(u)) < E(** Au *W A «* 

k=l 

and 

£([(4 v U *W] - i(jtk A u fc )JV£])A& j lj2(t k A u^iV^ALj 
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m 

+ 2 V u k )N%\ - [(t k A u fc )JVj^])A^. 

k=l 

Proof. Fix t, u > 0. Let v = (ti A ui, ti A u% . . . , t m A u m ) and w = 
(ti V ui , t2 V «2) • • • j tm V tt m ) . For l<i<m and 1 < J < I < N" , 

c ov(Xi,j(t),Xi,i(u)) = cov(l - X J l/ j(t), 1 - Xi,i(u)) 

with 

E[{i-xrj(t)Hi-xj; ( («)}] 

= n n E t(!- w 2 ] n 

fc=l ^ '=1 l={v k N%]+l > 

and nU(0 = lifo>6- 

Therefore after some straightforward algebra, we have that 

cov (xij(t)>Xi',{( u )) 

(3.4) =|n(i-^J K ^ i " [ ^ 1 } 

(m m \ 

n{( i - + Ku\ [vk ^ ] - ma - ^) a } [ ^ • 
fe=i fc=i / 

It follows trivially from (3.4) that cov(xf,-(t), Xiii u )) — 0- Since all the prod- 
ucts in (3.4) are less than or equal to 1, we have that 

cov( X r j (t),X^(u)) 

3.5 

m m 

fe=i fc=i 

Hence (3.2) is proved. 

Straightforward algebraic manipulation of (3.4) gives 

coy{Xi,j{t),xli(t)) - cov(x^(t),X^(u)) 

/ m m \ 

= mu - ^) 2 + a^}^-^ - n(i - M yK^i-M 

\fc=l fc=l / 

xcov( X ^(v),x^(v)) 

(3 - 6) 

+{m i -^) 2[ ^ ] } 
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x n{a-^) 2 +AL} [ ^ ] - K ^ 

\fc=l 

-n{(i-^) 2 > [ ^ ] - [ ^ ] ■ 

k=l / 

By arguments identical to those used in (3.5), 

m m 

Ilia - + Afc}^H"^ _ JI(i _ ^jMfl-MB 

fc=i 

m 

II{(l-^ a }^- [ ^-l 

=1 

m 

l_pj ( l_^.)K^]-«] 



fc=i 



(3.7) 



< 



k=l 
+ 



fc=l 



fc=i fe=i 

A similar inequality exists for the latter term on the right-hand side of (3.6) 
giving 



|cov(x l i ' 7 (t),xL(t))-cov(xr, 7 (t),xL(u))| 



(3.8) < 3 V u fc )JV£] - [(t fc A u fc )iV£])/4 N $>* A «*M<Mw< 

U=i J U=i j 

m 

fc=l 

Therefore, since the inequality in (3.8) holds with t and u reversed, (3.3) 
follows immediately. □ 

Lemma 3.2. For any s > 0, let D(s) = [0,s]. For 1 < % < m, and for any 
s> 0, 



(3.9) 

Thus 
(3.10) 



sup 
teD(s) 



1 



x»{t)- n {t) 



as v — > oo. 



sup |X 1/ (t)(^n)" 1 -r(t)| ^0 as oo. 
t6D(s) 
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Proof. Consider fixed 1 < i < m and t G D(s). Since rf (t) — >• j-j(t) and 
n% — ► 7r as v — ► oo, to prove that jspX" (t) — > ^(t), it is sufficient to show 

i 

that 



1 TX V 1 iVf 

— Xf (t) - ^-rf t) = -^-Xf(t - -*-rf (t 







as ^ — > ex). 



Fix e > 0. By Chebyshev's inequality, 



1 JVf , 



EE cov (^(t),xr/t)). 



Therefore, by (3.2) 



1 N u 

—xrw-^rm >e 



/c=i 



and so, for all t £ D(s) 



(3.11) 



J-xnt) - r<(t; 







as z/ — > oo. 



The remainder of the proof of (3.9) follows by arguments similar to [8], 
Lemma 5.1 since both X\{t) and n(t) are nondecreasing in t. Therefore the 
details are omitted. □ 



Corollary 3.3 follows immediately from (2.1), Lemma 3.2 and (3.1). 

Corollary 3.3. Suppose that for 1 < i < m, 7T U — > iv and C, u — > £ as 
v — ► oo. T/ien if £>0, 

1^ — >r as z^—>oo. 

Alternatively, if £ = 0, 

(3.12) min{|T^|,|T^-r|}^0 as z,^ oo, 

where t exists if and only if R > 1 . 

4. Branching process approximation. We shall assume that there exists 
a > such that for all v > 1, sl = a. (The extension to the case & u — > a as 
f — > oo is trivial.) Corollary 3.3, (3.12) is valid in this situation and therefore 
our aim is twofold: first, to find P(lim !/ _ +00 = r), that is, the probability 
of a major epidemic outbreak, and second, to find the limiting distribution 
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of conditional upon a small epidemic outbreak occurring (i.e., TJ^ 
as v — > oo). 

We shall show that converges in distribution to the total progeny Z 
of a suitably defined multitype Galton-Watson branching process. We pro- 
ceed by giving a slightly different construction of the epidemic process from 
that given in Section 2. It is straightforward to show that the two construc- 
tions are equivalent, with regard to the final size of the epidemic process. 
We shall then describe the branching processes to which the epidemic pro- 
cesses are coupled, along with the limiting branching process which has total 
progeny Z. 

We make the following assumption. For 1 < k < m, there exists a random 
variable XJ k such that 

V%N V IT U k = (U k>1 ,U k , 2) U Km ) as v - oo, 

with E[U kti ] = nia < oo (1 < i < m). For 1 < k < m, let R£(= (R kjl ,R kj2 , ■■, 

RIJ) = (Bin(iVr, V% A ),Bm(N%, V£ 2 ), Bin(JV£, F£J), with TV = (R?, 
R2,...,R„). Thus for 1 < i,k< m, let R ki denote the total number of 
initially susceptible individuals of type i contacted by a typical type-A; infec- 
tive. For 1 < i, k < m and I > 1, let PY k « { denote the set of type-i individuals 

contacted by the Ith. infective of type k, with R" k ; % i = \PK « J. Finally, let 
-pu r p^ p;/ pi/ \ 

Let Z^ = [Z\ , Z£ , . . . , Z^J denote the total progeny, excluding the a initial 
ancestors, of an m-type Galton-Watson branching process with (mixed bino- 
mial) reproductive law R^. That is, a typical type-A; individual has R ki off- 
spring of type i. For 1 < k < m, let R fc = (Po(U kA ), Po(U kj2 ), . . .,Po(U k>m )). 
Therefore let Z = (Zi,Z 2 , ■ ■ ■ , Z m ) denote the total progeny, excluding the 
a initial ancestors, of an m-type Galton-Watson branching process with 
(mixed Poisson) reproductive law, R = (Ri, R2, • • • , R m )- 

Lemma 4.1. Suppose that for l<k<m, V%N U W U fc as v — > 00 . 
Then: 

(i) R^ Rfc as v — > 00, 

(ii) T^-^Z as u^oo. 

Proof, (i) The result is a trivial extension of [13], Lemma 5.8(i). Hence 
the details are omitted. 

(ii) Note that an immediate consequence of (i) is that Ti v — ► Z as v — > 00. 
For < s G 7L m , let M v (s) denote the event that the infectives 

(l,l),(l,2),...,(l, Sl ),(2,l),(2,2),...,(2, S2 ),...,(m, Sm ) 
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make infectious contacts with distinct sets of individuals. The key point is 
that conditional upon M u (s) occurring, the epidemic and branching pro- 
cesses can be coupled, so that the processes coincide for the first s individ- 
uals. Now for s > 0, it is straightforward to show that 

(( m-1 m Si s k \ 
U U UU{ p M np (M)> 
I i=l k=i+lj=ll=l ) 

{m Si— 1 Si \ \ 

u u u ^/ p w)} H 
t=i j=i i=j+i ) i 

(4.1) 

m— 1 m Si Sfc 

<EE EEnpynp^^) 

i=l k=i+lj=l 1=1 
m Si-1 Si 

+ EE E P(PfonP^0). 

i=l j=l l=j+l 

For 7^ (k,l), 1 < d < m and 6, c > 0, 

6c 

P ( P (ij),d n P (k,l),d ¥= \ R (i,j),d = 6 ' R (k,l),d = C ) < ]m7 ■ 

Thus for / (fe,Z), 

P(P(ij) n P W) ^ ) ^ E Tn7E[l^ iJ))d ]E[flf fci0)d ] - 

Therefore the right-hand side of (4.1) converges to as v —> oo. 
Thus for all sG(Z+) m , 

P(T^<s) = P(T^<s|M^(s + a))P(iVr(s + a)) 

+ P(T^ < s\M u (s + a) c )F(M u (s + a) c ) 

= P(Z U <s\M u (s + a))¥(M u (s + a)) 

+ P(T^ < s\M u (s + a) c )P(M i/ (s + a) c ) 

^P(Z<s) asi/^oo, 

and the lemma is proved. □ 

By applying standard branching process theory (see, e.g., [15]), if R < 1 
(cf. Section 3), the total progeny of the Galton- Watson branching process 
is almost surely finite. In this case Lemma 4.1 is sufficient for studying 
T^,. Alternatively, if R > 1, there is a positive probability 1 — n*=i Qi' that 
the Galton-Watson process produces an infinite number of progeny, where 
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for s G [0, 1], hk(s) = E[njliexp((sj - l)itjU k j)] (l<k< m), and q is the 
unique solution in [0, 1) of q = h(q) = (/ii(q), ^(q), • • , h m (q)). 

Finally, using a lower-bound branching process approximation (see, e.g., 
[8, 9, 22]), it is straightforward to show that 

p (j™ T ~ = °) = n = 1 - p (jsa, T - = r ) • 

5. Central limit theorem for the final size distribution. In Section 3 we 
established a weak law of large numbers result for the final size distribution 
conditional upon there being a major outbreak (T^ — r as v — > oo). We 
shall extend the results of Section 3 to obtain a central limit theorem for 
(T^, — T)y/N u H as — > oo. We shall require the further constraints that for 
1 < hji k < m, \/ N u (N u [i"j — — > 0, and there exists < Ayfe < oo such 
that {N u ) 2 \"j k — > Xijk as z/ — ► oo. All standard epidemic models satisfy the 

above criteria (cf. Section 7). We also assume that \/N u (ir u — ir) — > as 
i/ — ► oo, with relaxations of this assumption considered in Section 6. 
For 1 < i < m and t > 0, let 

17(t) = yffi{jLx?(t) ~ f <(t)} 

and 

where P u = ITTr 1 and IF = diag«,vr£, . . . 

The general procedure in the literature for showing that (TJ^, — t)^N v H 
converges in distribution to a multivariate Gaussian distribution is as follows. 
First, show that the finite-dimensional distributions of Y v converge to the 
finite-dimensional distributions of a Gaussian process Y. The next step is to 
show that {Y^} are tight, and hence, Y u Y as v — > oo. Then by applying 

[11], Theorem 4.4, and pages 144-145, it follows that Y^(T^ + (?) 
Y(t + C) as v *■ 00 • Finally, it is shown that there exists an invertible 
mxm matrix J7 such that (T^ - t)VN u U and Y U (T^ + C")^" 1 have the 
same limiting distribution, namely, Y(r + ^)L r ~ 1 . 

However, we shall give a simpler and more direct proof that Y V (T% + 

(?) -^-> Y(r + £) as v — > oo. By [11], Theorem 4.1, it is sufficient to show that 

Y v (r + C) Y(t + C) and |Y^(T^ + C) - Y^(r + 01 as i/ oo. 
These results are proved in Lemmas 5.3 and 5.4, respectively. 

In proving a central limit theorem for {Y^t)}, we shall prove that Y^(t) 
has the same limiting distribution as Y i/ (t)P v +Y(t), where Y^(t) and Y(t) 
are independent and defined as follows. For 1 < k < m, and t > 0, let Z%(t) = 
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YT=\ Y}-=i n v$ J)tk , zi(t) = ET^iUN-m]^ and z k (t) = YT=i knmk, with 

Z"(t) = (ZZ(t),Z%(t),...,Z% l (t)), and z"(t) and z(t) denned in the obvi- 
ous fashion. For b € R m , let g(b) = (exp(— b\), exp(— 62), • • • ,exp(— 6 m )) and 
G(b) =diag(g(b)). Then 



(5.1) Y"(t) = -{g(Z"(t)) - g(z-(t))}v^n, 

and Y(t) is an m-dimensional Gaussian distribution with mean and covari- 
ance matrix (I — G(z(t)))G(z(t)). We shall therefore consider the limiting 
distribution of Y"(t) in Lemma 5.2, before considering Y u in Lemma 5.3. 
For t > and 1 < k < m, let 

_[tfe^] 

where /z^ = (/x^,//^, . . . ,/i^ m ). Therefore W£ fc describes the fluctuations 
about the mean infectivity of the first [ffciVT| infectives of type k. 

Lemma 5.1. Let Wt &e an m-dimensional Gaussian distribution with 
mean and covariance matrix V^n-fX^fcLi tk^k^-k}^/^, with = (Xkij) (1 < 
< m) . Then for t > 0, 

m m 

(5.2) W£ = ^{W^Vn} ^{Wt.feVn} = W t , say, as 1/ - 00. 

fc=l fe=l 

Proof. The proof is trivial, and hence the details are omitted. □ 

Lemma 5.2. For any t > 0, 

Y"(t) W t G(z(t)) as 00. 

Proof. By the mean value theorem, there exists = (s^, • • • > s m) 
lying between Z^(t) and z"(t) such that 

= -£ (ZJ(t) - ^(t))— g(y) ViVPIF. 

• =1 I y=s" J 



For 1 < j <m, 
d 

-j— g(y, 

d Vj 

where 6ij denotes the Kronecker delta. Thus 



= -((5i i exp(-Si),<52jexp(-4),.. . ,<5 mi exp(-s^)), 

y=s" 



Y"(t) = (Z"(t) - z u (t))G(s u )y/N u Tl 
= W£G(s"). 
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By Lemma 5.1, Z v (t) — z(t), and hence, s u — z(t) as v — > oo. Therefore, 
by [11], Theorem 5.1, Corollary 2, G(s") G(z(t)) as v — > oo. The lemma 
then follows from Lemma 5.1. □ 



Lemma 5.3. For any t > 0, 

Y"(t)-^Y(t) asu^oo, 

where Y(t) is an m- dimensional Gaussian distribution with mean and 
covariance matrix 

H(t) = {I-G(z(t))}G(z(t))+^^ 

Proof. Since P u —> I as v —> oo, Y(t) is the limiting distribution of 
Y u (t)P" + Y(t) as i/->oo. 

For 7 G R m , let Z*(t) = ££i 7/^(t), and for s e R, let 

A,(a;t;7)=E[exp(is^(t))]. 

Let 



/^(s;t;7) = J| exp 

k=l 

x E^ 



r fc (t){l-r fc (t)} 



J] exp( i S7fe ^K(t|V) - rUt)} 

k=l 



Therefore, using the Cramer-Wold device and characteristic functions, the 
lemma follows by proving that 

(5.3) |/„(s;t;7) - f u (s;t;j)\ asi/^oo. 

Conditional upon V 1 ", X^-(t) and x^,-/(t) are independent, and so, 

t; 7) = Ev„ [E z , [exp(M^(t))|V v ]] 



V" 





r N" 








exp( 




Lfc=ij=i 





exp ia 7fc _— =( x ^ ,(t) -rjfc(t)) 



For 1 < k < m, and 1 < j < N%, 



E-y" 



1 



exp is 7fe -=(xL(t) -rjfc(t)) 
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exp(* S ^(r£(t|V')-r£(t)) 



x E, 



expfw ™ {^(tJ-rfttlV")} 



V 



Let 



^(s;t; 7 )=Ev, 



Ill ex P 

fc=i 



7.S- 



fc ^ 7fc K(t|v-)-<(t)} 



x (i-£^i r -( t |v^){i-^(t|v^)} 

where vax(xfe,j| V 1 ') = r v k (t\V v ){l - r£(t|V")} (1 < k < m, 1 < j < JV£). 

For any s£l, and for any random variable W such that E[|VF| 3 ] < oo, 
there exists no £ N such that for all n > uq , 

E 



rap (, 5 -L^)]-(l + ^E W -l-E[^]) 



(5.4) 



Therefore for all 1 < < to, 1 < j < N% and t > 0, it follows from (5.4) that 



<(t=) n\wn 



|/ ! ,( S ;t;7)-/, 1 ( S ;t;7)|<E^f^) 3 -0 



fc=i 



as z/ — > oo. 



Let 



m r / ap 

exp * S ^7 fe K(t|V) - r£(t)} 



/ s 2 N v 

x exp(- y ^ 7 fr£(t|V''){l -r£(t|V)} 
It is trivial, using [12], page 94, Lemma 4.3, to show that 



|/>;t; 7 )-/>;t; 7 )l<I>£ 
For 1 < k < to, 



k=l 



2NI 



(5.5) 



r£(t|\>)- U-exp -E £ 

\ \ i=l j=l , 



as ^ — > oo. 



m [tiJVJ'] 

<EE 

i=i j=i 
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with E[ET=iE [ ff\v ( i Jlk n = JXiMTKAfo + (^) 2 ) - as oo. 

Therefore, by Lemma 5.1, and [11], Theorem 5.1, Corollary 2, 

(m 1*1^1 \ / m \ 

i=l j=l / V i=l / 

and so, r p k (t\\ v ) r k (t) as v -> oo. For all 1 < A; < m, iV£/iV£ -> 1 as 
— > oo and for all v > 1, |/£(s;t;7)|, |/„(s;t;7)| < 1. Therefore 

|^(s;t;7)-/„(s;t;7)| ->0 as z/ -> oo. 

Thus (5.3) is proved. □ 

For 1 < k < m, exp(—Zk(r + £)) = <7fc, and so, 

(5.6) G(z(t + 0) = S = diag(ai )C 7 2 , . . . , <r m ). 
Therefore, Y(r + C) ~ iV(0,H), where 

(5.7) H = {/ - £}£ + Ev^jf> fe + Cfc)vr fc A fc | yTIE. 

The next step is to show that Y v {T u OCJ + C) and Y u (t + Q have the same 
limiting distribution. 

Lemma 5.4. Suppose that -^U r. Then 

Y"(T£, + C) Y(r + C) as i/ - oo. 

Proof. By Lemma 5.3, Y"(t + C) Y(t + C) 

as ^ — > oo. Therefore, 
by [11], Theorem 4.1, to prove the lemma it is sufficient to show that 

|Y"(T£, + C) ~ Y"(t + C)|-^0 as v - oo. 

Fix e > 0. Then for 1 < i < m, by Chebyshev's inequality 

F(\Y?(T v 0O + C)-mT + Q\>e) 

= J2n\Y"(T V oo + O - KiT + C)| > e|T^ + C = k' / )P(T^ + C = k") 

k" 

< E 72 v ar(^(T^ + C) - Y?(t + C)|T£, + C = k")P(f + C = k") 
= lE[var(^(T^ + O - Y?(t + C)|T^ + C)}- 
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Thus 



P(|^(T^ + n-^(T + C)|> £ ) 

NY NY 

< -tt E E E[cov( X ^(T^ + O - xSj(t + 0, 



(5.8) 



X^(T^ + n-X^(r + C)|T^ + C 
r2 -{E[var( X r,i(T^ + O - xW + C)|TSo + 01 



+ (iVf - l)E[cov( X ^(t^ + C) ~ Xi,x(r + C), 

The lemma follows by proving that the right-hand side of (5.8) converges to 
as v — ► oo. 
First, 

E[var( X r,i(TSo + C) " X^(r + C)|T^ + C)] 
<E[|rr(T^ + C)-<(-r + C)|]- 

Since r, it follows that |rf (T£, + C") - rf(r + C)| as i/ -> oo. 

Therefore, since rl'(-) is bounded, the right-hand side of (5.9) converges to 
as v — > oo. 

By (3.3) of Lemma 3.1, 

NncovUl^ + O-xtiir + C), 

xWoo + C) ~ Xhij + C)|TSo + Ol 

(5-10) < 6iVr|E([(^, fe Vr fc M] - [(f^ k Ar k )m)^ 



.fe=i 

m 



+ 2N? £([CC, fe V t*)^ - A r fc )iV fc H)AL. 

fc=i 



It is trivial to show that for all 1 < k < m, 

[CC, fe V r fc )^] - [(T^ )fc A T k W k ]} -A as z, - oo. 



Therefore, since N v \i v ki — ► /i^ and (N u ) 2 X kii — > A^i as f — > oo, the right- 
hand side of (5.10) converges in probability to as z/ — > oo. Finally, since 
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0<T^<1, 



NfEilcovixl^ + O-xUT + C), 
as v — > oo. □ 



We are now in position to derive a Gaussian limit for (TJ^ — t)\/N u H. 

Theorem 5.5. Let U be the m x m matrix with entries Uij = 5ij — 
^TTiTTjfiijaj . Then U is invertible. Furthermore, if {C, u — £)V 'N u Ii v — > and 
— > t as v — > oo, then 

(T^, — r)\/iV^n —> N(0, (U T )~ 1 EU^ 1 ) as v -» oo, 
where S zs defined in (5.7). 

PROOF. Note that 

= (X^(T^ + Oi^Uy 1 - r(r + c))v^ri 



= Y"(T£, + C) + K(T^ + C)P V ~ r(t^ + Olv^TI 
+ {r(T^ + 0-r(T + C)}v / A^n. 

Hence, 



({T^ + C- r(T^ + C)} ~ W + C - r(r + C)})v^TI 

= Y V (T£, + C) + K(T^ + " r(T^ + 0}v^Tl 

+ {C - Clv 7 !^ 7 !!. 

By the mean value theorem, 

({T^ + C ~ r(T^ + C)} - {r + C - r(r + C)})v^TT 



t=s" , 



where s 1 ' lies between + and r + £. 

For t > 0, let <S(t) and £7(t) denote the m x m matrices with entries, 
Sij(t) =5ij-7TifJ,ij(l-rj(t)) and %(t) = — y / 7f77r7^ i j ( 1 — (t ) ) . Therefore 



CC - r)5(s") v^TI = (T^o - r) v^TlC/(s^ 

(5.11) 

= Y-(T^ + + A„ say, 
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where A u as v -> oo. Thus by [11], Theorem 4.1, (T^, - r)\/iV^n and 
Y iy (T^ + ^)L r (s 1/ )~ 1 have the same limiting distribution as v — > oo, provided 
that lim^oo C/(s^) _1 exists. 

Since ?/(■) is continuous, and s 1 ' — r + £ as ^ — > oo, it follows by [11], 
Theorem 5.1, Corollary 2, that 

[/(s 1 ') -^C/(t + C) = U asz/^oo. 

By [20], page 575, U is invertible. Therefore, by combining Lemma 5.4 and 
[11], Theorem 4.4, we have that 

y"(t^ + CMs'T 1 y(t + c)^" 1 ~ JV(o, (c/^^str 1 ) 

(5.12) 

as f — ► oo, 

and the theorem is proved. □ 

6. Random group allocation. There is an interesting variant of the above 
model in which individuals are independently and randomly allocated to a 
group according to a probability distribution, 7r. The main motivation for 
studying this variation of the model is that it enables us to consider a wider 
range of epidemics upon random graphs in Section 7. 

Let S denote an m-dimensional Gaussian random variable with mean 
and covariance matrix T = II — ir T ir. Then V N v (tv u — ir) — ► S as v — > oo. 

To assist in the analysis of the model, we define for 1 < % < m, X\ (t) = 

Y.% X^-(t) and fy = yffiijLjtttt) - r\ (t)), with A>f = N»m (cf. Section 

i 

3). Also X 1/ (t) and Y y (t) are defined in the obvious fashion. Note that N 1 ' 
is a random vector, while N 1 ' is deterministic. 

The results proved in Sections 3 and 4 are valid for the above model 
with only trivial modifications of the proofs required. However, Theorem 
5.5 requires minor alterations for the current model. We begin by consider- 
ing {Y u (T'^ JO + £")}. We then prove Theorem 6.2 which is an adaption of 
Theorem 5.5. 

Lemmas 5.1-5.4 hold with {Y V (T^ + £")} in place of {Y y (T^ + £")}. 
Therefore we shall require Lemma 6.1 to utilize these results for the current 
model. 

Lemma 6.1. 

|Y"(T£, + C) ~ Y^(T^ + 01-^0 asv - oo. 

Proof. Fix 1 < i < m and e > 0. Then by applying Chebyshev's in- 
equality in a similar manner to Lemma 5.4, we have that 

p(i^(T^+r)-^(T^+m> e ) 
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< -,E[var(l7(T^ + C) - 1?(T£, + OIN")] 



1 



' JV.f 



A > var E xij (TSo + c) - E (T^o + c 



£ 2 iV7 



-E 



E E ^(xUt^ + cuut^ + C)) 

■j=(N»AN»)+l l=(NVANy)+l 



Since TJ^ is bounded above by 1, it follows from (3.2) that 

F(\Y»(T» co + C)-Y»(T U oc + C)\>£) 
1 



(6.1) 



< 



< 



s 2 N? 
1 



-E 



\n?-n?\ i + |ivr-ivriE^ A 
i fc=i 



A « 



E 



7T; 



+ snp (AH 2 ALE 
Kfc<m 



« " "if 



7T; 



Thus the right-hand side of (6.1) converges to as v — ► oo, and the lemma 
is proved. □ 

Theorem 6.2. Suppose that (C - C)v / ^TT and T^, r as 
1/ — > oo . Then 

(T£, - r)v / iV^n iV(0, (E/ T ) -1 {3 + TjC/- 1 ) as i/ -» oo, 
where S and J7 are defined in (5.7) and Theorem 5.5, respectively, and 
T = v / lF T (/ - £)T (/ - S)vTI ZT . 

PROOF. The proof is similar to Theorem 5.5, and therefore only an 
outline of the proof is given. 

Following the proof of (5.11), there exists s u lying between T^, and r 
such that 

(T£, - t)Vn u UU(s v ) = Y^(T^ + C) + A, + (C - OV^H, 
where U(-) is defined in Theorem 5.5, and 

A u = {r"(T£, + C)P U - r(T^ + C)}VNni. 
Since r(t) = 1 - g(z(t)), let 



A„ = r(T^ + c)(ir - u)Vn»ii- 1 

= v^(7r" - 7r)(/ - G(z(T^ + C)))vTF T . 
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It is trivial to show that 

| A, -k u \ = |{r"(T£, + C) ~ r(T^ + C)}P u VNnFi | _L> (J as „^oo. 

Let A v = Vn^Itt" - tt}(J - tyVTF 1 . Since G(z(T^ + £")) G(z(r + 
£)) = S [see (5.6)] and \J~N^{t^ v — 7r} -^-> S as v — > oo, it follows that | A„ — 
Ajy| — > as — > oo. 

By Lemmas 6.1 and 5.4, and the triangle inequality, 

|Y"(T£, + C) - Y"(t + C)l as i/ ^ oo. 

Therefore 

(T^ - r)VWtt = {Y u (t + C) + A„)C/(s^)- 1 + M„, 

where Mj, — ► as oo. 

For all v > 1, Y^(r + £) and A„ are independent, with Y^(r + £) 
iV(0,H) and A„ S(/- E) vlF 1 as ^ — > oo. Therefore Theorem 6.2 follows 
since U (s u ) — ^ [/ as v — > oo. □ 

The adaption of Theorem 6.2 to other (random) group allocation schema 
is trivial. Theorem 5.5 is the special (deterministic) case where S = 0. 

7. Special cases. The generic results developed in the preceding sections 
can be specialized to a number of specific examples. These highlight the 
wide-ranging applicability of the model described in Section 2. 

7.1. Ball-Clancy (1993) model. This very general model is described in 
full detail in [6]. The important features from our perspective are the fol- 
lowing. There are m groups, and for v > 1, there are & u initial infectives 
and initial susceptibles. The susceptibles remain fixed within their ini- 
tial group. However, once infectious, an individual can move between the 
m groups. For 1 < i < m, let I* = I\, ■ ■ ■ , Z^,), where for 1 < k < in, 
P k denotes the time spent in group k by an infective who originates from 

group i, with P = J2T=i ^l- For 1 < hj, k < m, let denote the rate at 

which an infective, who originates from group i, and is currently in group j, 
makes infectious contacts with a given individual in group k. Let Bi = (f3 k j) 

(1 < % < m). Clearly for 1 < i,k < m, V» k = 1 - exp(-^ J2f=i Pkjfy, and 

hence, N£V» k » n k J2f=i £§4 Thus for 1 < i < m, 

V^jVTT Ui = VBfli as v -> oo. 

Therefore, if for 1 < i < m, E,[P] < oo and var(P) < oo (cf. [6], Section 4, 
page 727), it is straightforward to show that the conditions stated in Sections 
4 and 5 are satisfied. 
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7.2. Ball-Clancy (1995) model. In [7] an epidemic in a homogeneously 
mixing population with m different types of infectives is considered. An in- 
dividual on becoming infected chooses its type at random from {1,2,..., m} 
according to a probability distribution 7r, irrespective of the type of the in- 
fector. An equivalent construction is to randomly allocate a type to each 
individual according to 7v prior to the epidemic. Then should an individual 
become infected, it becomes an infective of its allocated type. Therefore this 
is a special case of the model described in Section 6. 

7.3. Random graph models. The main motivation for the current work 
is the study of the final size distribution of epidemics upon random graphs. 
A population (and epidemic) upon a random graph Q are defined as follows. 
Two individuals, i and j, are said to be acquaintances, if in the graph Q an 
edge exists between vertices i and j. 

A wide variety of extensions of the Bernoulli random graph can be shown 
to satisfy the randomized Reed-Frost criterion. We shall illustrate this range 
by considering static and dynamic multitype Bernoulli random graphs. The 
key property that we shall require is that the existence of an edge between 
any two vertices is independent of the remainder of the graph. This will 
allow the epidemic and random graph to be constructed in unison; see [16]. 

For a static multitype Bernoulli random graph Q, for 1 < i,j < m, let 
jjb (ctij > 0;N U > maxjajj}) denote the probability that there exists an 
edge between a given vertex of type i and a given vertex of type j. We 
shall require that the random graph is irreducible, that is, there is a nonzero 
probability of there existing a pathway of edges between any pair of vertices. 
Let Wij denote the probability that a type-i individual, while infectious, 
makes an infectious contact with a type-j acquaintance. Thus V£a = jftWij 
is the probability that a type-i infective, while infectious, makes an infectious 
contact with a type-j susceptible. Therefore the conditions of Sections 4 
and 5 are satisfied. The underlying multitype random graph can then be 
constructed in a similar manner to [16] with an edge existing between two 
vertices k and I, if in the epidemic individual k infects / or vice versa. An 
edge exists between the vertices k and I, according to the correct conditional 
distribution, if neither individual k nor I infects the other individual. 

An interesting special case of the static multitype Bernoulli random graph 
with random type allocation is the mixed Bernoulli random graph. The 
mixed Bernoulli random graph epidemic model assumes that the infectiv- 
ity, and connectivity, of all individuals are independent and identically dis- 
tributed according to random variables W and D, respectively, with W 
and D independent. That is, for individuals k and I with connectivity ran- 
dom variables Dj~ and Di, respectively, the probability that k and I are 
acquaintances is proportional to D^Di. Then given they are acquaintances, 
the probability that individual k, while infectious, makes infectious contacts 
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with individual I is Wk- We assume that there exist m£N, and 6, 7r 6 K m , 
such that J2iLi 7T i = 1 and P(D = 9j) = (1 < j < m). Therefore, indepen- 
dently at random assign individuals to one of m groups according to tv, that 
is, according to Di. Then for 1 < i,j < m, let an individual of type i and 
an individual of type j be acquaintances with probability 9i6j/N u . Thus 

Vft = OiOjW/N" and the conditions of Section 6 are satisfied. 

Epidemics upon random graphs with prescribed degree have received con- 
siderable attention in the literature; see, for example, [3, 4, 17]. That is where 
a vertex, i say, has Ei edges connecting the vertex to E{ distinct vertices. 
The random variables {E\, E2, ■ ■ ■ , E^v} are assumed to be (asymptotically) 
independent and identically distributed according to a nonnegative, integer- 
valued random variable E. Particular interest has been devoted to the so- 
called scale- free graphs [17], where the mean of E exists, but the variance 
does not. Unfortunately, such epidemic models do not permit a (multitype) 
randomized Reed-Frost representation since the existence of an edge be- 
tween any pair of vertices is no longer independent of the remainder of 
the graph. However, in most real-life scenarios it is realistic to assume that 
K[E] < var(E), and so, the mixed Bernoulli model with D chosen such that 
E[E] = K[D} 2 and var(E) = E[D] 2 (l + var(D)) will often provide a suitable 
alternative. (We have a scale-free mixed Bernoulli random graph if the mean 
of D exists, but the variance does not.) Thus it would be interesting to ex- 
tend the results for the mixed Bernoulli random graph to continuous and /or 
scale- free D. 

Finally, we consider dynamic Bernoulli random graph epidemics, which 
have previously been considered in [1] and [4]. For a random graph Q, and 
for 1 < i, j <m, consider two distinct vertices, one of type i, vertex k, say, 
and one of type j, vertex I, say. In the dynamic random graph, we alternate 
between there existing an edge between the vertices k and /, and there 
being no edge between the vertices k and I. The time from the formation 
of an edge between vertices k and until the dissolution of the edge is 
distributed according to S^j ~ Exp(p~ ). The time from the dissolution of the 
edge between vertices k and / until the reforming of the edge is distributed 
according to N u S^~j, where Sf^ ~ Exp(/9^). We shall assume that for 1 < 
i,j < pfj'Pij > 0, that is, edges can form and dissolve between all types 
of vertices. Thus for t > 0, the two individuals corresponding to vertices k 
and / are acquaintances at time t, if an edge exists between vertices k and 

/ at time t. Therefore, in equilibrium = + P li i/ _ is the probability that 

Pi] "r Pij 

the two individuals are acquainted at any given point in time. We shall 
assume that the dynamic Bernoulli random graph is in equilibrium. 

Only minor modifications are required for the following cases: no edges 
exist between certain types of vertices (all that we require is the graph Q 
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is again irreducible), the graph is part static and part dynamic and general 
probability distributions for S~j and Sfy 

For simplicity in exposition we shall consider a generalization of the model 
considered in [4]. For 1 < i < m, infected individuals of type i have indepen- 
dent and identically distributed lifetimes according to Qi, where E[Q?] < oo. 
While infectious (and acquainted), an individual of type i makes infectious 
contacts with an acquaintance of type j at the points of a homogeneous Pois- 
son point process with rate Pij. As noted in [4], the probability of more than 
one partnership between two individuals during the course of an infectious 
period is negligible. We therefore have the following two cases to consider: 
first, where the two individuals are acquaintances at the beginning of the 
infectious period, which we shall term initial acquaintanceships, and second, 
where the two individuals become acquainted during the infectious period, 
which we shall term secondary acquaintanceships. Suppose that an individ- 
ual of type i, and an individual of type j, are acquaintances at the point of 
time, t say, at which individual i becomes infected. Suppose that the indi- 
viduals remain acquainted until time t + Zij, where Zij = S^~j. Conditioning 
upon Qi = qi, the probability that an initial acquaintance of type j is infected 
is 1 — exp(— Pij(qi A Z^j)). The probability that during the infective's infec- 
tious period a link forms between the infective and an initially unacquainted 
individual of type j is 1 — exp(— j^pfjqi) ~ j^ptjQi- Conditional upon a link 
between the two individuals being formed, the time (relative to the start of 
the infectious period) at which the acquaintanceship is formed has probabil- 
ity density function g{x) = ^ exp(-xp+./A^)/{l - exp(-qipfj/N u )} « l/q { 
(0 < x < q^, and g(x) = (otherwise). Therefore the probability of infecting 
an individual of type j who is not an initial acquaintance is approximately 
1 m 

■J^JPijH J o V%{1 - exp{-/%((<?i - x) A S {j )}} dx. 
Hence, Vft = Af£ + o{l/N u ), where 

M?j = a£{l - eM-Pij(Qi A Z hj ))} 
Pij 



+ jfiL ^ ~ e M-Pij((Qi ~x)A Sr.))} 



Pij + Pij 



Therefore, for 1 < i < m, there exists a random vector U, with finite mean 
and covariance matrix such that V; 'N V W U,-. 
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